Ultrafast Spin Dynamics in Optically Excited Bulk GaAs at Low Temperatures 
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This paper presents a study of electron spin dynamics in bulk GaAs at low temperatures for 
elevated optical excitation conditions. Our time-resolved Faraday rotation measurements yield sub- 
nanosecond electron spin dephasing times over a wide range of n-doping concentrations in quanti- 
tative agreement with a microscopic treatment of electron spin dynamics. The calculation shows 
the occurrence and breakdown of motional narrowing for spin dephasing under elevated excitation 
conditions. We also find a peak of the spin dephasing time around a doping density for which, under 
lower excitation conditions, a metal-insulator transition occurs. However, the experimental results 
for high excitation can be explained without a metal-insulator transition. We therefore attribute 
the peak in spin-dephasing times to the influence of screening and scattering on the spin-dynamics 
of the excited electrons. 

PACS numbers: 72.25.Rb,71.70.Ej,78.47.-p,72.25.Fe 



I. INTRODUCTION 



Electron spin dynamics in semiconductors has experienced a dramatic revival of interest in the last decade due to 
the promise of semiconductor-based spintronics applications. 1 ' 2 In contrast to magneto-electronic systems involving 
metallic multilayers, semiconductors are better suited for the fabrication of integrated electronic and opto-electronic 
devices, but in non-ferromagnetic semiconductors the spin is not a conserved quantity due to the spin-orbit interaction. 
Thus an understanding of spin relaxation in semiconductors is essential, because it limits all information processing 
and storage using electronic spins. The renewed interest in combination with improved experimental techniques has 
led to important new results. For instance, long spin dephasing times in n-doped GaAs at low temperatures and 
low carrier densities have been measured using time-resolved Faraday rotation (FR) techniques. 3,4 Very recently, the 
technique of spin noise measurements has provided a new tool for the determination of "intrinsic" spin relaxation 
times without electrical or optical carrier injection. 5 These measurements show, for instance, the importance of a 
metal-insulator transition on the spin relaxation at low temperatures. However, the foundation for the theoretical 
explanation of these experiments above the transition temperature was already laid by an analysis of Dyakonov and 
Perel (DP) 6 who showed that in n-doped GaAs the decay of spin polarization occurs through the dephasing of the 
spin coherence in a spatially inhomogeneous effective magnetic field provided by the band structure. 

Although these "non-invasive" experimental techniques yield important results, one may also be interested in 
the excitation of appreciable spin-polarized carrier densities by optical or electronic means. Existing studies, which 
address these questions, have concentrated on the electron spin dynamics in n-doped GaAs-based semiconductors 
under electrical bias and Hanle effect measurements on a timescale of several hundred picoseconds. 8 

This paper is concerned with an investigation of the spin-dephasing dynamics in n-doped bulk GaAs after the 
excitation of spin-polarized carrier densities by ultrashort optical pulses. For a wide range of doping densities, we 
present a detailed comparison between time-resolved FR measurements and a microscopic theory, which is based on 
the original DP analysis. Our approach to compute the spin-dephasing time does not require the determination of 
auxiliary quantities, such as an effective momentum scattering time, which is often employed to achieve agreement with 
experiment. 9 Instead, the equations of motion for the spin and momentum-resolved electronic distribution functions 
are solved numerically including the relevant effects of the band-structure and the scattering of electrons with other 
electrons, phonons, and impurities. 10 Due to the numerical complexity of calculating the spin dynamics in bulk GaAs, 
theoretical results so far have concentrated on the microscopic spin dynamics in low-dimensional semiconductors, 11 ' 12 
with the exception of recent work, 13 or on accurate calculations of the analytical DP results for the spin decay. 14,15 

This paper is organized as follows. In Sec. II we explain the theoretical approach and present the dynamical equations 
for the spin-density matrix. In Sec. Ill we describe the experimental setup and our modeling of the experimental 
conditions in the framework of the theory. Sec. IV contains the experimental results and their interpretation by 
comparing them with model calculations. 
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II. THEORETICAL APPROACH 



We directly calculate the dynamics of the momentum resolved 2x2 spin-density matrix 

Pss ,{k) = (cl k c s , k ) (1) 



where ,s = f , j denotes the spin projection quantum number along the quantization axis z, c (c^) are electron destruction 
(creation) operators in the Heisenberg picture, and (• • • ) denotes the average with respect to the equilibrium statistical 
operator. The elements of the spin-density matrix are the carrier occupation numbers n s k = p ss (k) and the complex 
coherences = pn(k). In terms of the spin density-matrix, the average spin is given by 

(s x (k)) =^Re[* fc ] (2) 



Mk))=^(nl-ni). (4) 



The time evolution of p(k) is determined by the equation of motion 

(5) 
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where the coherent part includes the electronic energies, the influence of the effective magnetic field due to the band 
structure, and Hartree-Fock renormalizations. For the diagonal elements nf, we have 

J^felcoh = - j-{\^ B B + n x (fe)]im(* fe ) + n y (k)M*k)} 



andthc off-diagonal elements are given by 



^*feLoh = ^^bB + %n x {k) + n y (k)}(n{ - n k ) - % -n z {k)^ k 

+ 1 jX w *tt n l+* W»I - W 
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Here, s = ±1/2 denotes the electron spin and B = Be x in the Voigt geometry. The screened Coulomb interaction 

W q = ^— v q (8) 
in Eqs. (6) and (7) contains the bare Coulomb potential 

(9) 

and the dimcnsionlcss dielectric functions £bg and e q . The former describes the "background" screening due to the 
polarizability of the lattice, and the latter the screening due to the electrons provided by the doping. For this, we 
use the static limit of the Lindhard dielectric function 16 in the numerical evaluation. The scattering contributions 
(electron-electron, electron-phonon, and electron-impurity) to Eq. (5) can be obtained by equation-of-motion or 
Green's function techniques 17-19 following, e.g., Refs. 20-22. As an example we write out the expressions for Coulomb 
scattering and define the kernels for in- and out-scattering 

K k n l!q = (1 - <+,)< + ( X - <+<*)< - 2Rc(*^ 1+q * fel ) (10) 

and 

*C = <+,(! - <) + <+,(! - <) - 2Rc(*; i+q * fel ) (II) 
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so that the scattering contribution for the diagonal elements becomes 




scat 



^2 W q S ( e k - e \k+q\ + £|fei+q| - Efei) 



fcl.q 



x (K£ uq [(l - n%)n% +q - Re(*l* fc+ ,)] 



(12) 



For the off-diagonal elements we have 



fcl.q 




mass m = 0.067too- An important step for the numerical solution of Eq. (5) for bulk GaAs is to account for the 
strongly anisotropic intrinsic magnetic field ri(fc) by an expansion of p(k) in a set of orthonormal functions Li(d, <p). 
Such an expansion was already used in Ref. 23 but, in contrast to the original treatment, we choose an expansion in 
functions that describe the angular dependence of fl most accurately, namely: Li(d,tp) oc 1 and £2.3,4(1?, f) <x Q x ,y,z 
where the components of the intrinsic magnetic field are given by 



with a material specific constant 7. 24 

Before we turn to the numerical results from this theory, which describes nonequilibrium spin dynamics due to 
the Dyakonov-Perel mechanism, it should be noted that there exists a variety of effects which lead to relaxation and 
dephasing of spin polarizations in GaAs. However, in moderately n-doped GaAs at low temperatures under elevated 
optical excitation conditions, the Dyakonov-Perel mechanism is dominant. The quantitative agreement between the 
parameter free microscopic calculation and the Faraday rotation measurements additionally, presented in Sec. IV, 
supports this point of view. For definiteness, we mention the following competing processes: Besides the intrinsic 
magnetic field the g- factor dispersion of electrons in GaAs 2 leads to additional dephasing. We neglect this effect 
because the optically excited spin polarization is due to electrons in a small energy interval above the Fermi energy 
(less than lOmeV). In this region the spread of rotation frequencies due to the g- factor dispersion Agfj,s\B\ is typically 
one order of magnitude smaller than the internal magnetic fields h{\fl(k)\). Thus the dephasing due to the g-factor 
dispersion is of minor importance for the theoretical analysis of the experiments presented in this paper. The Bir- 
Aronov-Pikus mechanisms usually acts on a longer timescale for the relatively low density of optically excited holes 
and therefore does not play a dominant role, either. 25,26 Relaxation due to interactions with a nuclear spin polarization 
is not included in the calculation; in the experimental results it is excluded by periodically switching the hclicity of 
the circularly polarized excitation beam. 



We have performed ultrafast FR measurements to determine the electron spin dynamics in n-doped bulk GaAs 
layers. 1,27,28 The doped GaAs wafer was mounted on a sapphire substrate and thinned down to about 30 ^m for 
the transmission measurement. Spin polarized electrons were optically excited by a circularly polarized pump pulse. 
In this setup, the ultrafast laser pulses are generated by a mode-locked Tksapphire laser oscillator operating at 
76 MHz repetition frequency. The samples were held at low temperatures in an optical cryostat with a split-coil 
superconducting magnet, and a magnetic field of typically 4T was applied. After the excitation, the electron spins 
begin to precess in a magnetic field applied parallel to the sample surface (Voigt geometry). By measuring the 
polarization rotation of a linearly polarized probe beam of the same wavelength transmitted through the sample, the 
electron spin dynamics can be monitored directly. The spin precession results in an oscillating signal, as shown in 
Fig. 1, that decays exponentially. We use an exponential fit for the signal envelope to extract the decay time due to 
dephasing processes, which is conventionally denoted 2 by the ensemble spin dephasing time T 2 *. 



n y (k) 
n z (fc) 



jk x (k 2 y -k 2 z ) 

n fk y (k 2 — fc 2 ) 
lk z (k 2 x - fc 2 ) 



(14) 
(15) 
(16) 



III. EXPERIMENT 



4 




200 400 600 

Delay Time (ps) 



FIG. 1: Time-resolved Faraday rotation signal for Nd op = 3.6 x 10 16 cm" 
(top), and numerical calculation (bottom). For the calculation iV exc = 1 
T 2 * is determined by the exponential decay time of the envelope. 



3 , B = 4T, and T — 4K as obtained by experiment 
x 10 16 cm -3 was assumed. The spin-dephasing time 



For the comparison with the experiment, we use for the calculation initial carrier distributions n = n^° v + n| xc , 
which contain both the influence of the doping and the optical excitation, because the excitation pulse duration is 
much shorter than the timescale of the subsequent electronic spin dynamics, so that any coherences introduced by 
the pulse can be ignored for the electronic dynamics. The itinerant carriers introduced by the doping are assumed 
to be completely relaxed, and arc modeld by unpolarized Fermi-Dirac electron distributions n dop = f(e%) at lattice 
temperature with electronic density equal to the density of dopants. We model the electronic carrier distribution n| xc 
excited by the ultrashort circularly polarized laser pulse propagating in z direction as a Gaussian with spectral width of 
the pulse, a = 15 meV. In accordance with the experimental procedure, we take the photon energy fko p approximately 
equal to the electron-hole transition energy at the Fermi momentum of the doping density under consideration. As 
in the experiment, we include a magnetic field of 4T, a lattice temperature 4K. Further we choose a 50% initial spin 
polarization of the optically excited carriers, which corresponds to the optimal spin polarization achievable by optical 
orientation at the bandgap. 23 For the calculation of the time-resolved FR signal shown in Fig. 1, we assume a linearly 
polarized probe pulse degenerate with the circularly polarized pump. Assuming that the holes are unpolarized 29 ' 30 
on the timescale of interest and that excitonic effects do not contribute, the FR signal is determined by 31 

e F cx]r ff (fc)[4-4] (i7) 

k 

Here, g(k) oc exp[(e£ — e\ — Fuj p ) 2 /<r 2 ], with a the spectral width of the pulse and e\ = fi 2 fc 2 /(2mh), mi, = 0.457mo, 
contains the contribution of the electron- hole transition at wavevector k to the signal at the probe photon energy ftw p . 
For the experimental and theoretical "raw data," as shown, for instance, in Fig. 1, we use the same exponential fit 
for the decay of the signal envelope to extract T 2 * . 



IV. RESULTS 



Figure 2 shows the dependence of the spin-dephasing time T 2 * on the photoexcited carrier density. Both in theory 
and experiment we find monotonically decreasing dephasing times for n-doped GaAs (A*d op = 3.6 x 10 16 cm~ 3 ) and 
increasing dephasing times for undoped GaAs with the better agreement for the undoped GaAs. When comparing 
measured and calculated results, it should be noted that there is a considerable uncertainty in the experiment because 
the photoexcited carrier density is difficult to estimate, and for the doped sample, we have an additional uncertainty 
in the doping density. Taking this into account, the agreement for the doped GaAs is very good. The agreement for 
undoped GaAs is even better because the weak dependence of the spin-dephasing time on the excited carrier density 
means that the experimental uncertainty in the optically excited carrier density has only a limited influence. 

The spin-dephasing times for undoped and n-doped GaAs in Fig. 2 differ by more than a factor of three and have 
different slopes due to the different spin dynamics in both cases: In the undoped material, the spin of a noncquilibrium, 
highly spin-polarized (50%) electron plasma relaxes, whereas in the doped material the spin dynamics of polarized 
optically excited electrons occurs in the presence of unpolarized electrons residing in the material due to the n-doping. 
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FIG. 2: Measured (symbols) and calculated (lines) spin dephasing times vs. photoexcited carrier density for doped (upper 
curves) and undoped (lower curves) bulk GaAs at T = 4K and with B = 4T external magnetic field. The doping density is 
3.6 x 10 16 cm" 3 (calculation) and (3.6 ± 0.3) x 10 16 cm" 3 (experiment). 

For these two scenarios, the interplay of the different scattering mechanisms (carrier-carrier, carrier-phonon, and 
carrier-impurity) and screening, and thus the influence of the optically excited carrier density on the spin-dephasing 
time are different. 

Figure 2 already shows the pronounced influence of the doping density on the spin-dephasing time, and we investigate 
it now in more detail. The theoretical and experimental results in Fig. 3 are obtained by varying the doping density 
for a fixed density of optically excited carriers of 2 x 10 16 cm -3 . We find good agreement of the measured data points 
with our microscopic calculation including only itinerant-electron dynamics, which suggests that the peak of the spin- 
dephasing time at intermediate doping densities can be explained without invoking a metal-insulator transition. 7 A 
qualitative explanation of this peak is provided by the motional narrowing relation 2 ' 32 



In Eq. (18), (f2) is the effective Larmor frequency due to the DP spin splitting averaged over the electrons taking 
part in the spin-precession and dephasing dynamics. The spin and momentum relaxation times are denoted by T sp i n 
and t p , respectively. For low and moderate doping densities, i.e., for densities up to 10 17 cm -3 we are in the motional 
narrowing regime because the Larmor frequency fi(fc) of the precession of an electron is proportional to fc 3 . Thus 
it becomes large only for electrons excited at larger momenta, which is the case at higher doping concentrations 
for optical excitation beyond the quasi Fermi-level of the relaxed electrons. The combination of elastic impurity 
scattering and mainly inelastic Coulomb scattering leads to a decrease of the effective momentum scattering time r p 
with increasing doping density. Via the motional narrowing relation (18) this translates to a longer spin-dephasing 
time. For higher doping concentrations, the motional narrowing regime described by the inequality in Eq. (18) is left 
and the spin-dephasing time drops with increasing doping concentration. 2 

For a more quantitative understanding of the doping density dependence in Fig. 3, we analyze the contributions of 
the different scattering mechanisms. To this end, we also show in Fig. 3 the spin-dephasing time for the two artificial 
cases of (i) only elastic impurity scattering, and (ii) no impurity scattering. For case (ii) we still include carrier-carrier 
and carrier-phonon scattering, which lead to clastic and inelastic scattering events, but the inelastic processes have 
the larger available phase space and therefore dominate case (ii). For undoped GaAs, i.e., doping densities of less than 
10 15 cm" 3 , the full calculation is identical to case (ii), which shows that carrier-carrier and carrier-phonon scattering 
are responsible for the spin dephasing. 35 For intermediate doping densities between 10 15 cm -3 and 10 17 cm" 3 the rise 
of the spin-dephasing time towards a maximum at 8 x 10 16 cm" 3 is determined by the interplay of the elastic impurity 
scattering, the predominantly inelastic carrier-carrier/phonon scattering, and the influence of screening on all three 
scattering processes. Inelastic or elastic scattering processes alone cannot explain the behavior of the spin-dephasing 
time in this region: If only elastic scattering is included, as in case (i), electrons are not scattered away from the 
higher k states, in which they are excited by the laser; there they experience a faster precession, cf. Eq. (14), and thus 
dephase more quickly with increasing doping concentration. This effect counteracts the increasing effectiveness of the 
carrier-impurity scattering, which, in the picture of Eq. (18), leads to a shorter r p . If only carrier-carrier and carrier- 
phonon scattering is included, as in case (ii), the optically excited carrier density of 2 x 10 16 cm" 3 dominates over or 



Tspin C*- 



— if (ft)Tp < 1. 
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FIG. 3: Spin dephasing time vs. doping density at 4 K and with 4 T external magnetic field. Shown together with the calculated 
result including all scattering mechanisms (solid line) are curves obtained by including only impurity scattering (dashed grey 
line) and without impurity scattering (dash-dotted grey line) . The triangles are experimental values. The excited carrier density 
for calculation and experiment is always iVexc ~ 2 x 10 16 cm -3 . 

is comparable to the doped carrier density for doping concentration of up to several 10 16 cm~ 3 , so that no pronounced 
density dependence is to be expected for these scattering mechanisms alone. 36 For high doping concentrations beyond 
10 17 cm -3 , the optically excited carriers cannot be scattered out of the high k states in which they are created, because 
of the high density of relaxed electrons. Then, by Eq. (14), the precession of the excited electrons becomes too fast for 
the timescales of elastic [case(i)] scattering, inelastic scattering [case(ii)] or their combination [full result], so that the 
spin-dcphasing time experiences a pronounced drop. 2 The complex behavior of the doping-density dependence of the 
spin-dephasing time is thus quantitatively explained by the interplay of the different elastic and inelastic scattering- 
mechanisms, and the influence of screening. This picture is essentially supported by a reanalysis of the data of Fig. 3 
in Ref. 33. Last, but not least, we wish to stress that the numerical results provide a microscopic picture of motional 
narrowing in spin dephasing, i.e., an increasing T 2 * with decreasing momentum scattering time, and the breakdown 
of this behavior for fast electronic precession. We do not need to resort to established analytical results for DP spin 
dephasing, 23 which do not work over the whole density regime and for pronounced non-equilibrium conditions. 34 

V. CONCLUSION 

We presented results for the spin dynamics of itinerant electrons due to the DP spin relaxation mechanism for a 
wide range of doping concentrations, and under elevated excitation conditions. Our calculations yield quantitative 
results for the DP spin relaxation process, which show the occurrence and breakdown of motional narrowing in 
spin dephasing for elevated optical excitation conditions. We find a maximum of the spin dephasing-time on a sub- 
nanosecond timescale at intermediate doping concentrations close to the density where, at low excitation conditions, 
a metal- insulator transition occurs. By numerically investigating the different contributions to the spin-dcphasing 
time under elevated excitation conditions, we showed that the maximum of spin dephasing times in the intermediate 
doping regime is due to the interplay of elastic and inelastic scattering mechanisms, as well as screening. 
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